Nuclear spin dynamics in the quantum regime of a single-molecule magnet 
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We show that the nuclear spin dynamics in the single-molecule magnet Mni2-ac below 1 K is 
governed by quantum tunneling fluctuations of the cluster spins, combined with intercluster nuclear 
spin diffusion. We also obtain the first experimental proof that - surprisingly - even deep in the 
quantum regime the nuclear spins remain in good thermal contact with the lattice phonons. We pro- 
pose a simple model for how T-independent tunneling fluctuations can relax the nuclear polarization 
to the lattice, that may serve as a framework for more sophisticated theories. 
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Single-molecule magnets (SMMs) are nanometer-sized 
high-spin molecular clusters organized in a crystalline ar- 
ray, which sets the direction for the anisotropy axis 0. 
Reversal of the cluster spin can occur either classically, 
by thermal activation, or quantum mechanically, by tun- 
neling through the barrier Q, Q ■ SMMs are attractive 
model systems to study the effects of coupling magnetic 
qubits to the environment (nuclear moments, phonons), 
with the associated problems of decoherence and the lim- 
its of quantum mechanics at the large scale 0, • For 
both aspects the hyperfine coupling between cluster spin 
and nearby nuclear spins is expected to play a crucial but 
subtle role: since this coupling is many orders larger than 
the quantum tunneling splitting, a static hyperfine in- 
teraction completely blocks tunneling. Contrariwise, by 
considering it as a dynamic bias that sweeps the electron 
spin levels through the tunneling resonance, Prokof'ev 
and Stamp (PS) have argued that this interaction in fact 
promotes incoherent tunneling events 0. 

Experimentally, although time-dependent magnetiza- 
tion experiments [7j showed a \fi dependence and isotope 
effects agreeing with the PS predictions [6j, fundamental 
aspects of the spin-dynamics, like the essential role of nu- 
clear spin diffusion, remain to be verified. Further, in the 
PS model the quantum relaxation of the cluster spin is to 
the nuclear spin bath and is expected to be many orders 
of magnitude faster than conventional spin-lattice relax- 
ation to phonons. A crucial test is thus whether or not 
the experimental nuclear polarization relaxes to the lat- 
tice (phonon) temperature, even at such low T that only 
electron spin tunneling fluctuations are left ("quantum 
regime" ) . Interestingly, whereas in order to relax to the 
lattice the nuclear spins generally need electron spin fluc- 
tuations, in the quantum regime those same nuclei would 
provide the only source for such fluctuations via the PS 
nuclear-spin-mcdiatcd quantum tunneling model. It is 
by no means obvious how such a T-independent process 
could establish thermal equilibrium between spins and 
lattice. Here we report a NMR study of the dynamics of 



55 Mn nuclei in [Mni 2 Oi2(02CMe)i 6 (H 2 0)4] (Mni 2 -ac), 
which experimentally answers the above points and poses 
a crucial test for a realistic description of the coupling be- 
tween a magnetic qubit and its environment. 

Mni2-ac is the SMM with the highest anisotropy bar- 
rier (~ 65 K) discovered so far; its core is composed 
of 4 Mn 4+ ions (electron spin s = 3/2), and 8 Mn 3+ 
ions (s = 2) in two inequivalent crystallographic sites. 
The intracluster superexchange interactions lead to a to- 
tal spin S — 10 for the cluster. Below T w 3 K the 
electron spins are effectively frozen along the anisotropy 
axis, thereby enabling 55 Mn NMR even in zero applied 
field, by exploiting the local hyperfine field -Bhyp felt by 
the nuclei. This allows the use of nuclear spins as local 
probes for the fluctuations of the cluster spin by study- 
ing the nuclear spin-lattice relaxation (NSLR) and the 
transverse spin-spin relaxation (TSSR), without disturb- 
ing the zero-field tunneling resonance. We have chosen 
the resonance line of the 55 Mn nuclei in Mn 4+ ions, hav- 
ing a central Larmor frequency uj^jl-K sa 230 MHz and 
a relatively small quadrupolar splitting The exper- 
iments were performed on Mni 2 -ac crystallites, cast in 
Stycast 1266 epoxy and oriented in 9.4 T magnetic field 
at room temperature. The NMR coil with the sam- 
ple was placed inside the elongated tail of the plastic 
mixing chamber of a specially designed dilution refrig- 
erator. This allows a continuous flow of 3 He around 
the sample and assures excellent thermalization. The 
NMR signal was detected by spin-echo technique, with 
typical duration t^/ 2 = 10 fj,s for the 90° pulse. Since 
the 55 Mn nuclei have spin I — 5/2, the recovery of the 
nuclear magnetization, M z (t), after an inversion pulse 
obeys 0: M z {t)/M z {oo) = 1 - [(100/63) exp(-30Wf) + 
(16/45) exp(-12Wi) + (2/35) exp(-2Wi)], where W is 
the NSLR rate [13 [Fig. 1(c), solid lines]. The TSSR rate 
T 2 _1 is obtained by a single exponential fit of the decay of 
transverse magnetization, M xy (t) = M xy (0) exp(— i/T 2 ), 
except at the lowest T where also a gaussian com- 
ponent T^q needs to be included, yielding M xy {t) = 
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FIG. 1: (a) T-dependence of the NSLR (•) and TSSR (■) 
rates in zero field and ZFC sample, (b) Decay of echo inten- 
sity at T = 20 mK in ZFC (▼) and FC (V) sample. Inset: 
in the ZFC sample half of the nuclear spins (black arrows) 
has Larmor frequency —uin instead of +ujn because of the 
reversed orientation of the cluster spin (gray), (c) Recovery 
of nuclear magnetization after an inversion pulse, at T = 20 
mK in ZFC (a) and FC (A) sample. 



A/4 y (0)exp(-t/T 2 )exp[-0.5(t/T 2G ) 2 ] [Fig. 1(b), solid 
lines] . 

Between 1 and 2 K, both the NSLR and the TSSR 
show a roughly exponential T-dependence [dashed curves 
in Fig. 1 (a)] , which is well understood in terms of the 
fluctuations of -E>h yp produced by thermal activation of 
the electron spin levels [HEEl. The NSLR rate can 
be obtained from the spectral density at u = ujn of the 
transverse component of the fluctuating part of -Bhyp, 
with the implicit assumption that -Bhyp fluctuates around 
its average direction (which coincides with the molecule's 
anisotropy axis) but does not flip over, as in a tunneling 
event. Extrapolating the observed high-T NSLR to the 
mK range would lead to astronomically long relaxation 
times. In a preliminary work |l3| we observed that, upon 
cooling down to 20 mK, the NSLR saturates to a roughly 
T-independent plateau, indicating that only fluctuations 
due to quantum tunneling within the ground doublet are 
contributing to the relaxation. The crossover between the 
thermally activated regime and the quantum regime [Fig. 



1(a), solid lines] is clearly visible at 0.8 K, in agreement 
with magnetization experiments [LH . 

The value W « 0.03 s" 1 of the NSLR found below 
0.8 K is surprisingly high, considering that the relax- 
ation of the global magnetization in Mni2-ac takes years 
at low T. On the other hand, it is well known that 
any real sample of Mni2-ac contains a fraction of fast- 
relaxing molecules (FRMs), which are characterized by 
one or two distorted local anisotropy axes for the Mn 3+ 
ions 0]; for those molecules the barrier is reduced to 
35 K or even 15 K [ltij . yielding much faster tunneling 
dynamics. At the same time, however, we have verified 
that the observed NMR signal comes from nuclei in stan- 
dard, slow-relaxing molecules, even though the electron 
spin of such molecules remains frozen during the experi- 
ment. The fluctuating dipolar field produced by a tunnel- 
ing FRM on the nuclei of neighboring (frozen) molecules 
is far too small to account for the observed NSLR, so 
we have suggested |L3J that the relaxation mechanism 
should involve intercluster nuclear spin diffusion (not in- 
cluded in Refs. 0,0]), linking nuclei in frozen molecules 
to those in FRMs. By studying the magnetization depen- 
dence of T 2 _ 1 we can now provide strong evidence for the 
proposed mechanism. When comparing the TSSR in a 
demagnetized, zero-field cooled (ZFC) sample (where the 
cluster spins are randomly oriented up or down) with a 
saturated, field-cooled (FC) sample (where all spins have 
the same direction), we find that the FC sample has a 
faster TSSR, with a ratio T 2G 1 (FC)/T 2G 1 (ZFC) « 1.35, 
very close to V2 [Fig. 1(b)]. In terms of intercluster spin 
diffusion this has a simple explanation: in a FC sample 
the nuclei in equivalent crystallographic sites of differ- 
ent molecules have the same Larmor frequency, lun, thus 
flip-flop transitions are possible with all neighbors. In a 
ZFC sample the nuclei are divided in two groups having 
Larmor frequencies +u>n or — u>n, depending on the local 
spin orientation; for the nuclear dipole-dipole interaction 
this is equivalent to having diluted the FC system by a 
factor 2, yielding a \[2 times smaller TSSR [l7| [see inset 
of Fig. 1(b)]. The presence of a predominantly gaussian 
component as found in the TSSR at low T confirms the 
importance of nuclear dipolar couplings. 

Further insight in the relationship between the dynam- 
ics of the central quantum spin and the nuclei is provided 
by the field dependence of the NSLR. Applying an exter- 
nal field B z parallel to the anisotropy axis destroys the 
resonance condition for tunneling, thereby hindering the 
fluctuations needed for the NSLR; this explains the pro- 
nounced peak in W{B Z ) found around zero field as shown 
in Fig. 2(a). In comparing ZFC and FC sample, it is seen 
that both the width of the resonance and the zero-field 
value are quite different. In particular, one may con- 
clude that there are more tunneling events at zero field 
in the ZFC sample, as could be seen already from the 
difference in nuclear inversion recovery [Fig. 1(c)]. Such 
an observation, which is obviously impossible to obtain 
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FIG. 2: (Color online) (a) Longitudinal field dependence 
of the NSLR rate W{B Z ) at T = 20 mK in the ZFC (•) 
and FC (o) sample. The measuring frequency was set to 
lu n (B z )/2tt = 230 + 10.57B, MHz. (b) W(B Z ) in ZFC sample 
at T = 20 mK (T) and T = 720 mK (a): a small peak is visible 
at the first level crossing B z as 0.5 T only at the highest tem- 
perature. In this dataset we used u>n{B z )/2tt = 231 + 10. 57-B z 
MHz, which better matches the center of the NMR line at high 
T. The insets show a sketch of the electronic level scheme with 
the observed transitions. 
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FIG. 3: Comparison of bath temperature Tbath (solid lines) 
and nuclear spin temperature T nuc i (circles), while cooling 
down (a) and while applying step-like heat loads (b). The 
waiting time between NMR pulses (see inset) was 60 s in (a) 
and 180 s in (b). Both datasets are at zero field in ZFC 
sample. 



by means of "macroscopic" magnetization measurements, 
should provide a critical test for more detailed models of 
the NSLR. The signature of tunneling fluctuations at the 
first level crossing around B z « 0.5 T, i.e. when the spin 
states with S z = +10, +9, . . . come in resonance with 
S z = —9,-8,... , becomes visible as a small peak in 
W(B Z ) only upon warming up to T — 0.72 K, i.e. close 
to the border with the thermally activated regime [Fig. 
2(b)] 0. 

Finally, we address another essential aspect of the 
dynamics of the coupled system of nuclear and cluster 
spins, so far not studied theoretically or experimentally, 
namely: "what is the nuclear spin temperature"? In 
other words, is the nuclear spin polarization indeed re- 
laxing to an equilibrium value dictated by the lattice 
phonons, which are in thermal contact with the 3 He bath 
at temperature Tbath? In that case the intensity of the 
NMR signal as a function of temperature should obey 
the Curie law M Z (T) = K/T. The calibration factor 
K can be defined at the highest T by assuming that 
there the nuclear spin temperature T nuc i equals Tbath, 
and then be used to convert the NMR signal intensity 
into an equivalent T nuc i while cooling down the system. 
As shown in Fig. 3(a), we find that T nuc i indeed follows 
the time evolution of Tbath, the small discrepancy below 
0.2 K being most probably due to heating effects of the 
NMR pulses. Data taken with a lower pulse rate [Fig. 
3(b)] demonstrate that, even below 0.1 K, the nuclei al- 
ways closely follow the evolution of Tbath- This direct 
experimental proof of an energetic contact between nu- 
clear spins and phonons confirms earlier results from the 
(field-dependent) low-T specific heat |l9|. in which siz- 



able amounts of nuclear and electron magnetic entropy 
were observed to be removed below 0.5 K. Since there 
is no relevant direct energetic coupling between phonons 
and nuclei, the thermalization of the nuclear spin system 
must involve the interplay with the electron spins and 
their coupling to the lattice, even in the T-independent 
quantum tunneling regime. 

A basic question to answer is what happens to those 
nuclei that belong to a molecule where a tunneling event 
takes place (in our case a FRM) , assuming that the neigh- 
boring molecules are frozen. For the ease of discussion, 
we shall consider N nuclear spins 1 = 1/2 per cluster, 
subject to a hyperfine field -Bhyp parallel to the anisotropy 
axis of the molecule: the latter assumption simulates the 
real situation for 55 Mn in Mni2-ac. The standard way 
of calculating the rate of transition between nuclear Zee- 
man levels as a consequence of a perturbing fluctuating 
field is useless here, since the Zeeman levels themselves 
completely change after each electron spin flip, so per- 
turbation theory is not applicable. A more realistic ap- 
proach is to recall that each electron spin level is split by 
hyperfine interactions into a quasi-continuum manifold 
of levels 0, H3 that can be labeled by the local nuclear 
polarization AN = — N*, which yields a hyperfine 
bias £jv (typically ~ 0.1 K). Since the hyperfine fields 
before and after the flip of the cluster spin are just an- 
tiparallel, the manifolds of Zeeman levels on cither side 
of the anisotropy barrier are simply the mirror of each 
other. Moreover, since the tunneling traversal time is 
much shorter than l/u>jsr, the probability that a nuclear 
spin would coflip with the electron spin is negligible. This 
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implies that the only relevant tunneling transitions are 
those that do not require any nuclear coflip, thus AN 
= const. Considering the small additional bias £d 
due to dipolar fields from neighboring cluster spins, the 
tunneling transition with AN = const, requires an initial 
hyperfine bias such that = £,d- Once the molecule has 
tunneled, the hyperfine bias becomes £/v — — £d since the 
nuclear polarization is unchanged but -Bhyp is reversed; 
the new local hyperfine energy can then be redistributed 
to other nuclei via intercluster spin diffusion until the 
equilibrium within the nuclear spin bath is achieved. In 
this way the effect of tunneling is the "conversion" of 
dipolar into hyperfine energy and vice versa. 

Our data show that this description is still insufficient: 
to obtain a nuclear magnetization in thermal equilibrium 
with the lattice, tunneling events must be accompanied 
by creation/annihilation of phonons. In our opinion a 
crucial role may be played by the Waller mechanism 
[2l| . i.e. the change in the dipolar field when the dis- 
tance between neighboring molecules is modulated by 
lattice vibrations. Even at very low T we can expect 
the existence of low-energy phonon modes that corre- 
spond to displacements of the clusters with respect to 
each other. Here we consider the cluster cores as rigid 
objects within the soft matrix (with Debye temperature 
Qd ~ 20 K) of the ligand molecules, an approach success- 
fully used to account for the Mossbauer recoil-free frac- 
tions of metal cluster molecules • The modulation of 
the dipolar field, whereby the total bias may sweep back 
and forth through the tunneling resonance, can thus pro- 
vide a probability of incoherent tunneling with emission 
or absorption of phonons, whose energy would be released 
or extracted from the nuclear spins in the way described 
above, i.e. using the tunneling of electron spins as in- 
termediary. The detailed balance between emission and 
absorption may then provide the equilibration of nuclear 
spin and lattice temperatures. 

In conclusion, we have shown that the nuclear spin dy- 
namics in Mni2-ac below 0.8 K is driven by tunneling 
fluctuations of the cluster electron spin, in combination 
with intercluster nuclear spin diffusion and thermal equi- 
librium between nuclear spins and phonon bath; the lat- 
ter aspect calls for the extension of existing theories of 
incoherent quantum tunneling within the ground doublet 
to include inelastic processes. 
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